BMC Molecular and Cell Biology

BackgroundLamins are major regulators of the spatial and functional organization of chromatin. Lamins at the nuclear periphery form the lamina that anchors heterochromatin to the nuclear envelope. A subpool of A-type lamins localizes in the nuclear interior, where they also bind to euchromatic genomic regions. A-type lamin properties and chromatin association are regulated by lamin-associated polypeptide 2alpha (LAP2). Here we systematically analyze, how LAP2 depletion affects chromatin organization, accessibility and gene expression on a genome-wide level. ResultsLAP2 depletion in mouse dermal fibroblasts positively and negatively affects chromatin accessibility and gene expression throughout the genome, which correlates with changes in chromatin association of A-type lamins and the nucleosomal remodeler proteins BRG1 and CHD4. In particular, A-type lamins bind to open chromatin regions close to BRG1 and CHD4 binding sites and deregulated genes, but do not directly accumulate on genes and BRG1 and CHD4-enriched sites. Unsupervised clustering of the datasets on LAP2-bound genomic regions confirms spreading of A-type lamins to active chromatin regions containing deregulated genes and an enrichment of chromatin remodelers on a subset of these genomic regions. ConclusionsLAP2 depletion in fibroblasts leads to a gross rearrangement of chromatin. Genome-wide chromatin reorganization is linked to spreading of A-type lamins to active chromatin regions and accompanied by a restriction of chromatin remodelers to a subset of active genomic regions. These changes correlate with changes in chromatin accessibility and gene expression throughout the genome, particularly in regions where lamin binding is gained in LAP2 knockout versus wildtype cells.

Mutations in mitochondrial DNA (mtDNA) compromise ETC activity and impair oxidative phosphorylation. Since eukaryotic cells contain multiple copies of mtDNA, the resulting phenotype depends on the proportion of mutant mitochondrial genomes (the heteroplasmy level). Using isogenic cell lines carrying similar mtDNA deletions, a linear decline in cellular respiration was observed as mitochondrial DNA heteroplasmy increased. Despite this, cellular redox imbalance did not change until heteroplasmy exceeded 50%. As heteroplasmy increased past 70%, cells also exhibited an integrated stress response (ISR) and impaired translation was observed. These defects were reversed by either addition of asparagine or overexpression of pyruvate carboxylase (PC). The dependence on exogenous asparagine in other respiration-deficient cells was found to correlate inversely with the PC expression level. For example, patient-derived thyroid tumor cells, harboring high heteroplasmy for a Complex I mtDNA mutation and low levels of PC, exhibited asparagine auxotrophy, and L-asparaginase treatment suppressed tumor growth. Together, these findings demonstrate a role for mitochondrial pyruvate carboxylase in cellular asparagine synthesis under conditions of compromised respiratory activity.

The methyltransferase EZH2 (Enhancer of Zest Homolog 2) and the demethylase KDM6B (Lysine Demethylase 6B) have been associated with epithelial to mesenchymal transition (EMT) and poor prognosis in various cancers. These enzymes methylate and demethylate H3K27me3 and regulate distinct sets of genes controlling EMT induction, despite having opposite catalytic activities. This could be due to their recruitment or the modulation of their activity by partner proteins on specific loci. This work sought to identify proteins associated with chromatin and interacting with EZH2 or with KDM6B during EMT. To do so, co-immunoprecipitation and mass spectroscopy was used under TGF{beta} (Tumor growth factor {beta}) and TNF (Tumor necrosis factor ) treatment to induce EMT in A549 lung cancer cells. Surprisingly, numerous proteins related to focal adhesions were identified to interact with EZH2 or KDM6B. These proteins are part of a nuclear protein interaction network previously described as nucleo-adhesome. Among these proteins, TGFB1I1 (transforming growth factor induced peptide 1) and CSRP2 (cysteine and glycine rich protein 2) were further confirmed to interact with KDM6B in the nucleus and even more so during EMT. The target genes of these complexes were then sought by knocking down KDM6B, TGFB1I1 or CSRP2. Three genes (coding Integrin alpha 5, Laminin y2 and Matrix Metalloproteinase 9) were confirmed to be regulated by KDM6B, TGFB1I1 and CSRP2. These findings may have clinical relevance, as immunohistochemistry analyses performed on a cohort of lung cancer patients revealed increased nuclear localization of TGFB1I1 and CSRP2 in cells undergoing EMT.

Injuries to individual cells happen frequently as a result of physiological and environmental stresses during their normal daily functions that can lead to a ruptured cell cortex (plasma membrane and underlying cortical cytoskeleton). The capacity of cells to rapidly repair general daily injuries, as well as ones resulting from trauma, infection, or diseases/cancer, is essential for their survival. While we know the general cell biological outline of the highly-conserved physiological events taking place during cell wound repair, our knowledge of the molecular mechanisms governing the repair process is still fairly limited, due in large part to the lack of information regarding the molecules, machineries, and pathways involved. Here, we performed a genetic screen of 1322 fluorescent-tagged proteins to identify cell wound repair components that are recruited upon laser wounding or whose expression is lost and/or altered upon laser wounding. We identified 129 proteins that are recruited to wounds during the cell repair process through high resolution spatio-temporal expression analyses of these gene fusions in conjunction with a fluorescent actin reporter. Strikingly, we find that many members of the Rab family GTPases are recruited to wounds where, in addition to their well-known roles in intracellular membrane trafficking, they are affecting actin cytoskeletal organization and dynamics during the repair process. These studies are allowing us to define the earliest acting proteins, as well as those required at specific steps in the repair process based on their recruitment patterns and the precise timing of their recruitment to wounds. Thus, our imaging-based screen is providing us with a global view of the repair processes, as well as a large number of genes/gene families that provide new entry points for examining specific steps in the cell wound repair process. Author SummaryCells in our bodies get injured every day from normal activity, environmental stress, infection, or disease. To survive, they must quickly repair these injuries and restore normal function. While some molecules have been identified as key players of cell wound repair, many of the molecules involved and their roles remain unknown. In this study, we identified new molecules that are involved in different steps of cell wound repair. Using laser-induced injury in the Drosophila model, we examined 1322 proteins and observed their spatial and temporal dynamics in a cell after injury. From the 1322 proteins examined, we identified 129 proteins recruited to distinct regions around the damage site during cell wound repair, suggesting roles in specific steps of the repair process. Interestingly, a subset of these proteins are Rab family GTPase members, highlighting new roles for these proteins in regulating actin dynamics. By identifying new candidate repair molecules, we provide a foundation for understanding how cells maintain their integrity and how repair processes may be influenced by factors such as wound size, infection, aging, and disease.

Gasdermins (GSDMs) are a family of pore forming protein that trigger pyroptosis by permeabilizing cell membranes. Pyroptotic cells often release the proinflammatory cytokines interleukin-1{beta} (IL-1{beta}), and IL-18, thereby promoting an inflammatory response. GSDMs are typically cleaved by caspases or granzymes, which enable their translocation to the membrane. Here, we showed GSDMA and GSMDB localize to the cytoskeletal fraction of keratinocytes. Disruption of cell junctions causes gasdermin A and B (GSDMA and GSDMB) to translocate to the membrane fraction in the absence of cleavage. Cell junction disrupted keratinocytes release post-translationally modified keratins, but not IL-1{beta} or IL-18. These events depend on endocytic mechanisms associated with recycling of cell junctional proteins. Our study suggests that cell junction disruption can drive translocation of GSDMA and GSDMB from cytoskeleton to plasma membrane in keratinocytes, however there may be a subsequent trigger that causes the confirmational change allowing these gasdermins to form open pores.

Endothelial-to-mesenchymal transition (EndMT) has become a central mechanism in developmental biology, fibrosis, vascular disease, and cancer. We previously reported on an integrated signaling model in which TGF-{beta}2 induces EndMT through coordinated activation of Smad-dependent and Smad-independent signaling pathways converging on Snail, while GSK-3{beta} regulates Snail activity. We performed a systematic figure-by-figure reproducibility analysis of the original publication. Independent studies published between 2011 and 2026 were identified and curated according to predefined inclusion criteria. Each original experimental conclusion was evaluated for independent confirmation. In parallel, selected biochemical experiments were independently reproduced using newly acquired reagents and contemporary Western blot methodologies. Independent publications consistently reproduced each major mechanistic conclusion of the original study, including activation of Smad, ERK, PI3K/AKT, and p38 MAPK signaling, regulation of Snail expression, EndMT-associated marker switching, and GSK-3{beta}-dependent control of Snail activity. Independent laboratory experiments reproduced the principal biochemical findings using contemporary reagents and experimental workflows. The combined literature analysis and independent laboratory replication demonstrate that the mechanistic framework in our previous study has remained reproducible across multiple laboratories, endothelial cell types, disease models, and fifteen years of investigation. This work illustrates a complementary framework for assessing reproducibility that integrates direct experimental replication with cumulative independent validation.

RhoG is a member of the Rho-family of small GTPases, and is closely related to the canonical Rac1 GTPase, implicated in membrane trafficking, dorsal ruffling, macropinocytosis, and cell protrusion, but its activity has been difficult to visualize directly in living cells with high spatial and temporal resolution. Here, we developed and validated a genetically encoded, single-chain Forster resonance energy transfer (FRET) biosensor for RhoG based on a C-terminal full-length RhoG and an intramolecular RhoG-binding domain derived from ELMO1. The biosensor showed a robust dynamic range when comparing constitutively active and inactive RhoG mutants, responded appropriately to regulation by RhoGDI, GAPs, and GEFs, and detected growth factor-stimulated RhoG activation in live cells. Imaging in mouse embryonic fibroblasts revealed dynamic RhoG activation at leading-edge protrusions, dorsal ruffles, and forming pinocytic and macropinocytic structures. To define the signaling relationship between RhoG and its closely related family member Rac1, we combined the RhoG biosensor with a near-infrared Rac1 FRET biosensor for simultaneous live-cell imaging. Morphodynamic mapping showed that both RhoG and Rac1 activities were positively coupled to edge protrusion, with strongest correlations near the leading-edge, but their direct coupling varied with distance from the edge, indicating partial spatial decoupling within protrusive regions. Inhibition of Src-family kinases altered RhoG dynamics, strongly suppressed Rac1 coupling to protrusion, and inverted the normal positive correlation between RhoG and Rac1 activities. Signaling microdomain analysis further showed that Src inhibition selectively prolonged Rac1 microdomain lifetimes without significantly affecting RhoG domains. Together, these results establish a new biosensor for direct visualization of RhoG activity and reveal that RhoG and Rac1 are coordinated but spatially and temporally distinct components of protrusion-associated signaling networks, with Src-family kinases playing a central role in maintaining their normal coupling.

Immune checkpoint blockade has revolutionized cancer therapy, but the therapeutic efficacy is limited. Clinical trials on blockade of newly identified immune checkpoints didn't show promising result, suggesting that it might be insufficient to understand the function of immune checkpoints in cancer merely in the context of immunity. Here, we found mutually exclusive expression patterns of the immune checkpoint VISTA (or VSIR) and the neural stemness factor SETDB1, an oncoprotein that promotes immunoevasion, in xenograft tumors, suggesting that cells with high VISTA expression represents a differentiated, and hence, less or non-malignant state in tumor. Non-neural differentiation factors HHEX, MYOD1 and PPARG promote, whereas oncoproteins KRAS (and the mutant KRAS(G12D)) and SOX2, both being embryonic neural factors, repress VISTA expression. This tendency can be inferred from the finding that neural stemness is the core property of cancer cell. Manipulated expression of VISTA in cancer cells generated no significant effect on cell tumorigenicity and differentiation state, but led to change in cell morphology and actin cytoskeleton. Mechanistically, VISTA regulates a key cytoskeleton regulator, WASF2, leading to the change in cell morphology, which might interfere with signal transduction of immune response. The results suggest that 1) high expression of a protein in tumor might represent a less or non-malignant state, targeting of which would leave malignant cells intact, and consequently, leading to weak or even no therapeutic efficacy, a key factor worth considering for target selection; 2) immune checkpoints might play other roles in cells that interfere with regulation of anti-tumor immunity.

Compared with skin wounds, oral mucosal wounds heal more quickly, with minimal scarring, faster re-epithelialization, and reduced inflammation. One differentiating factor may be the differential transcription factor-associated gene networks involved in tissue regeneration. One such transcription factor, BATF3, was recently shown by us to promote wound-healing responses in vitro and in vivo. Our prior analyses also suggest that CREB5 is a differentially regulated transcription factor in oral wounds and may be involved in early wound-healing gene expression programs. CREB5 expression was induced in immortalized skin keratinocytes (HaCaT) to examine its effect on in vitro wound healing relative to immortalized gingival keratinocytes (TIGK). CREB5 overexpression let to differential expression of predicted downstream genes and improved skin keratinocyte migration in vitro. This work suggests that examining transcription factors and gene networks that regulate wound-healing responses in the oral mucosa may lead to the discovery of novel targets to improve skin wound healing.

Increased temperature is one of the first environmental cues encountered by bacteria upon entering a mammalian host. Here, we investigated the effects of temperature on the transcriptome and proteome of Escherichia marmotae and E. coli. Previous studies demonstrated that temperature affects motility in E. marmotae; therefore, we examined how temperature alters gene expression at 37 {degrees}C versus 28 {degrees}C and whether this response is conserved in E. coli. Strains were grown under static conditions at both temperatures, and gene expression and protein abundance were assessed by RNA transcriptome analysis and global proteomics. Temperature altered the expression of 111 genes (2.7%) in E. marmotae and 99 genes (2.5%) in E. coli (adjusted p < 0.05, [&ge;]2-fold change), with changes concentrated within specific functional pathways. In E. marmotae, flagellar and chemotaxis genes and operons involved in cellulose-dependent biofilm formation and nitrate respiration were markedly downregulated at 37 {degrees}C. In contrast, genes associated with fimbrial adhesion and immune evasion, including fimA/fimB, ompT, and prophage-associated loci, were upregulated. Proteomic analysis corroborated these trends, showing reduced flagellar and chemotaxis proteins and increased stress-adaptation and host-interaction proteins. E. coli showed a distinct response, with stronger enrichment of metabolic and amino-acid biosynthesis pathways and minimal changes in motility regulation. Together, these findings demonstrate that E. marmotae motility is temperature-dependent and may represent a mechanism for immune evasion within the host.

In TMPRSS2 fusion-positive prostate cancer, ERR is involved in regulation of ERG and promotes the androgen receptor independent signaling in the cancer progression. The ERR binds to the ERREs (estrogen-related receptor response elements) present at -5042 bp of the TMPRSS2- promoter and enhances the ERG overexpression that causes prostate cancer progression. To dissect the structural basis of the ERR recognition to the TMPRSS2 promoter DNA, we have purified the full-length ERR (ERRFL), NTD deleted construct (ERR{Delta}NTD), and the DNA-binding domain (ERRDBD) proteins and performed the binding analysis with 30 bp TMPRSS2-promoter DNA (5' -AGTCCAAGGTCGGTGGATC ACAAGGTCAGG-3'). Circular dichroism analysis showed that all three ERR proteins adopt native secondary structures. DNA binding induced subtle changes in the secondary structures, while enhancing the thermal stability (Tm) of all ERRa proteins. Binding analysis showed that ERRDBD bound weakly to the DNA, whereas ERRFL and ERR{Delta}NTD exhibited substantially higher affinities ~120-fold and ~131-fold than ERRaDBD, respectively. Small-angle X-ray scattering (SAXS) analyses revealed a dimeric ERRFL structure and an ERRFL-DNA complex (2:1) structure in solution and fitted well with Alpha Fold model of apo and DNA bound complex of ERRFL. Furthermore, 100 ns dynamics simulations on apo and DNA-bound ERRa proteins showed that all proteins remained structurally stable, with flexibility largely confined to loop regions of ERRa proteins. Our biophysical, DNA binding and structural analyses have revealed the mechanism involved in ERR recognition of the TMPRSS2- promoter DNA, which provides insight into ERR-mediated transcriptional regulation and development of anticancer drugs against ERR-driven prostate cancer.

Mutations in p53 and its 12 isoforms can alter its functions. As N-terminally truncated isoforms of p53 (delta40p53, delta133p53, and delta160p53) participate in tetramer formation, they are important regulators of cancer fate. Although delta40p53- and delta133p53-mediated regulation of cancer is well reported, the mechanism underlying delta160p53 production and its functional role remains unclear. We investigated the internal ribosomal entry site (IRES)-mediated translation of {Delta}160p53 and its role in cancer regulation. As differential synthesis of delta160p53 was observed under different stress conditions, IRES-mediated translation of this isoform was demonstrated using bicistronic luciferase constructs. No cryptic promoters or splicing sites were detected in the IRES sequence. Cell death and late apoptosis were significantly decreased, while proliferation, the number of cells in the S phase, and drug resistance were induced by delta160p53. Furthermore, delta160p53 did not induce p53-responsive promoters. RNA sequencing analysis of delta160p53 overexpression showed similar results, along with the inhibition of other tumor suppressor genes. Overall, our results provide insights into IRES-mediated translation of delta160p53, which can be considered a novel target for cancer treatment.

The mechanical properties of the nucleus are critical for maintaining nuclear integrity and function. We previously showed that chromatin dominates short-extension mechanics whereas lamins provide long-extension strain stiffening. To distinguish the roles of lamin isoforms, micromanipulation nucleus force measurements were performed on isolated nuclei from lamin A/C (Lmna-/-) and lamin B1 (Lmnb1-/-) knockout mouse embryonic fibroblast cells. Lamin A/C knockout does not alter short-extension nuclear stiffness but is essential for strain stiffening at longer extensions. Oppositely, lamin B1 loss reduced short-extension stiffness due to facultative heterochromatin loss while long-extension strain stiffening was slightly increased. Loss of lamin A/C and B1 resulted in similar lamin-chromatin linkers effects as LBR did not change and LAP2{beta} decreased in both. A simulation model of a polymeric lamin shell with stiff lamin A/C and softer lamin B1 subunits can qualitatively recapitulate experimental measurements of lamin knockout cells. Lamin A/C knockout resulted in abnormal nuclear shape but not nuclear blebbing or rupture whereas lamin B1 knockout, similar to other perturbations that cause heterochromatin loss, resulted in increased nuclear blebbing and rupture. This work illuminates the distinct mechanical roles of lamin A/C and B1 in determining nuclear structure and integrity.

THAP9 is a transposable element-derived gene which encodes a protein that is homologous to the active Drosophila P-element transposase (DmTNP). Both THAP9 and DmTNP possess a C-terminal domain (CTD) which is functionally uncharacterized. Sequence and structural analysis suggest that the THAP9-CTD has a novel fold which is only found in THAP9 homologs. To explore the evolutionary history and characteristics of this novel domain, exhaustive phylogenetic analysis (using MSA, structure prediction, MSTA-based clustering) was performed. THAP9-CTD homologs were more widely distributed throughout the animal kingdom in comparison to DmTNP-CTD homologs which were restricted to arthropods. Moreover, the THAP9-CTD homologs were more conserved, especially among mammals and birds and their average length increased in a class-specific manner. Comparison with the DmTNP-CTD homologs demonstrates that although their respective CTDs may have evolved independently, they both surprisingly share similar secondary structure elements consisting of three conserved helical regions made of hydrophobic residues that are predicted to make up a conserved core. The role of the respective CTDs were further investigated by creating truncation mutants lacking the CTD. Interestingly both THAP9 and DmTNP truncation mutants are still capable of DNA excision and integration suggesting that their respective CTDs are not essential for DNA transposition. Moreover, CTD truncation favours DNA integration in THAP9: this suggests that CTD acquisition during evolution may have led to THAP9 domestication as observed in other transposable element-derived genes like Rag1 and piggybac, which have similar terminal regulatory domains.

Gene model for the ortholog of Density regulated protein (DENR) in the Apr. 2013 (BCM-HGSC Dpse_3.0/DpseGB3) Genome Assembly (GenBank Accession: GCA_000001765.2) of Drosophila pseudoobscura. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Complex rearrangements are one of the rarest types of structural variants (SVs) and can be divided into two categories: complex chromosomal rearrangements (CCRs) and complex genomic rearrangements (CGRs). CCRs include structural rearrangements that present at least three breakpoints and show exchange of genetic material between more than two chromosomes and CGRs are rearrangements that present more than one junction and/or more than one SV in cis. They are usually formed by one of the chromoanagenesis mechanisms, where a massive disruptive cellular event leads to multiple structural rearrangements. Classical cytogenomic techniques have been commonly applied for their characterization, but methodologies that involve longer DNA molecules, namely optical genome mapping (OGM) and long-read genome sequencing (lrGS), present a considerably higher SV detection resolution, revealing more details about the rearrangements, including precise breakpoint location. Here, we describe six patients with complex rearrangements investigated through a combination of different techniques: karyotyping, chromosomal microarray, and OGM were performed to characterize the rearrangements. Subsequently, lrGS was used to further resolve the alterations, refine their breakpoints' location, and sequence their junction points. Three patients presented CCRs involving three, four, and six chromosomes, while three exhibited CGRs involving one different chromosome each, providing a variety of complex SVs to show the importance of each technique and their combination in rearrangement resolution. In total, the complex rearrangements presented 127 breakpoints, 66 junction points and involved 14 of the 24 chromosomes. Higher-resolution techniques revealed additional complexity in all cases. Despite the advances provided by OGM and lrGS, conventional karyotyping remained indispensable for complete rearrangement resolution. In two patients, the findings supported a novel mechanism combining features of the different chromoanagenesis processes. Furthermore, evidence of inherited alterations was identified, and the comprehensive characterization of the rearrangements enabled more accurate genotype-phenotype correlations. Our findings indicate that an integrated approach combining karyotyping, OGM, and lrGS can completely resolve SVs, including complex rearrangements.

Colorectal cancer cells increase glycolysis to help meet the metabolic demands required for cell growth. Many factors, both endogenous and exogenous, likely drive cellular metabolism and enhance glycolytic flux in colorectal cells. Interleukin-1 beta (IL-1{beta}) is a pro-inflammatory cytokine that is elevated in colorectal cancer. In this study, we investigated the effect of IL-1{beta} toward driving the cancer cell to increase glycolysis, while also suppressing the oxidation of the fiber-derived nutrient butyrate. The results presented here demonstrate that IL-1{beta} stimulated glycolysis and inhibited maximal mitochondrial respiration. IL-1{beta} also increased the phosphorylation of AKT and hypoxia-inducible factor 1 alpha (HIF1) levels. Utilizing colorectal cancer cells with AKT1/2 or HIF1 knocked out showed the requirement of these proteins in mediating the increase in glycolysis following IL-1{beta} treatment. Importantly, AKT1/2 was identified as upstream of HIF1, as IL-1{beta} still increased phosphorylation of AKT even in the absence of HIF1. However, loss of AKT1/2 completely abolished the ability of IL-1{beta} to increase HIF1 protein levels. Tumor necrosis factor alpha (TNF), another cytokine found to be elevated in colorectal cancer, also increased glycolysis in an AKT and HIF1-dependent manner. Our data point to a common pathway through AKT activation and HIF1 upregulation, by which pro-inflammatory cytokines increase glycolysis in colorectal cancer cells to help promote cancer progression.

MLLT1 (also named ENL) is a chromatin reader protein whose encoding gene was originally identified as a chromosomal translocation partner with MLL(KMT2A) in acute leukemia. However, its role in normal hematopoiesis has not been investigated. This study uncovers a critical role of Mllt1 in normal B cell lymphopoiesis. We found Mllt1 to be essential for early B lymphocyte development using a conditional Mllt1 knockout mouse model that we developed. A significant decrease of bone marrow B-lineage progenitors, splenic transitional B cells and peripheral blood B cells were observed in Mllt1del mice compared to control Mllt1fl/fl mice. Similarly, Mllt1 deletion in in vitro cultured B-enriched progenitor cells from Mllt1fl/fl; Rosa26CreERT2/+ mice resulted in reduced B cells, demonstrating the cell-intrinsic role of Mllt1 in this process. Direct MLLT1 target genes including Il7r and critical B-lineage transcription factors, Ebf1 and Pax5, were decreased following Mllt1 deletion. Gene set enrichment, gene ontology, and functional analyses of Mllt1-deficient cells showed significant alterations related to B cell development, critical relevant signaling pathways, DNA replication, and mitochondrial function. In vitro complementation with MLLT1 rescued the B cell phenotype observed with endogenous Mllt1 deletion; however, specific MLLT1 YEATS domain mutants lacking chromatin reader and RNA-binding functions were unable to rescue the phenotype. Taken together, our research demonstrates a previously unappreciated role for MLLT1 as critical for maintenance of B cell lymphopoiesis.

The human bony labyrinth has attracted increasing interest because of its taxonomic, evolutionary, and functional significance. Although sexual dimorphism has been reported in several aspects of the temporal bone, the extent to which sex, age, size, and allometry contribute to labyrinth shape variation remains poorly understood. Here, we investigated patterns of sexual shape variation in the human bony labyrinth using three-dimensional geometric morphometrics in a sample of 98 archaeological individuals from South America with known genomic sex. Centroid size and allometric effects were assessed in a subset of 90 individuals with comparable metric scaling. In addition to analysing the complete labyrinth, the cochlea and semicircular canals were examined separately to evaluate region-specific patterns of sexual shape variation. Principal Component Analysis showed extensive overlap between females and males, and overall labyrinth shape did not differ significantly between sexes. Males exhibited significantly larger labyrinths than females, and centroid size explained a small but significant proportion of overall shape variation. Regional analyses showed no evidence of significant sexual shape differences in the cochlea or in any individual semicircular canal when analysed separately. In contrast, the combined semicircular canal system exhibited subtle but significant sexual shape variation independent of centroid size, whereas morphological disparity did not differ between sexes. The geometric comparison of the female and male consensus configurations further showed that sexual shape variation was regionally heterogeneous. Whereas the cochlea exhibited a pattern of localized changes with low directional coherence, the semicircular canals displayed more coordinated regional shape changes. The male consensus also exhibited slightly higher canal circularity across all three semicircular canals, particularly the posterior canal, while differences in canal-plane orientation remained minimal. These findings demonstrate that sexual shape variation in the human bony labyrinth is subtle and anatomically partitioned among its components. Although significant sex differences in centroid size were detected across most anatomical regions, overall labyrinth shape and cochlear morphology were primarily influenced by allometry, whereas significant sex-related shape differences were detected only when the semicircular canals were considered as an integrated anatomical system. These findings demonstrate that sexual dimorphism in the human bony labyrinth is subtle but regionally heterogeneous, with the cochlea and semicircular canals exhibiting distinct patterns of shape variation, suggesting that these structures are influenced by different developmental, functional, and evolutionary processes.

Gene model for the ortholog of Insulin-like peptide 3 (Ilp3) in the D. pseudoobscura Apr. 2013 (BCM-HGSC Dpse_3.0/DpseGB3) Genome Assembly (GenBank Accession: GCA_000001765.2) of Drosophila pseudoobscura. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.